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Abstract. 

The propagator of the discrete Schrodinger equation is computed and its properties 
are revealed through a Feynman path summation in discrete space. Initial data 
problems such as diffraction in discrete space and continuous time are studied 
analytically by the application of the new propagator. In the second part of this 
paper, the analogy between time propagation and 2D scattering by ID obstacles is 
explored. New results are given in the context of diffraction by edges within a periodic 
medium. A connection with tight-binding arrays and photonic crystals is indicated. 
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1. Introduction 

The formulation of Quantum Mechanics in discrete variables is relevant to various 
aspects of undulatory physics. The applications range from solid state (tight-binding 
arrays p], polymers [2], photonic structures [3]) to fundamental questions concerning 
the structure of space (discrete canonical variables [I], cellular automata [5], algebraic 
deformations and bound states [6]). The dynamics of this kind of systems seems to 
be a natural area of study, filling a gap between the energy domain description and a 
time domain formulation yet to be explored. Moreover, among the many dynamical 
problems that one may pose, we distinguish diffraction by edges and the propagation 
of discontinuities (in space and in time [7]) as the most prominent representatives. In 
continuous space, these examples contain the effects of infinitesimal wavelengths. In 
discrete space, the natural frequency cut-off might be used to conjecture new effects. 
But beyond these conjectures, it is necessary to solve the propagation problem of specific 
distributions. 

In this paper we tackle some time- dependent problems by computing explicitly 
the propagator of an infinite tight-binding chain. This is equivalent to the discretized 
one-dimensional space propagator or the Green's function of a two-dimensional photonic 
structure in the lowest frequency band. The mathematical tools employed here, although 
of a simple nature (Bessel functions), constitute an original contribution of this work. 
They shall be thouroughly developed in order to get new results in the problem of 
diffraction. This applies to propagation in a periodic background and the effects 
produced by a minimal spacing between points. 

The paper is structured as follows: In section 2 we compute the propagator for the 
discrete Schrodinger equation by spectral decomposition (one or many dimensions for 
cubic lattices). As this object differs from its continuous counterpart, it is necessary to 
reveal its properties and its connection with Feynman paths in space and time. This 
shall be done in sections 2.1 and 2.2 respectively. In section 3 we move to specific time- 
dependent problems by applying the propagator to distributions of the edge type: A 
point-like distribution viewed as a minimal length wavepacket (section 3.1), the discrete 
version of the Moshinsky shutter (section 3.2) and a poissonian profile generalizing the 
former examples (section 3.3). Section 4 establishes the connection between diffraction 
in time in discrete space and the scattering of two-dimensional waves by a screen. Here 
it is shown that the scattering problem in a periodic background (such as that provided 
by a photonic structure) can be viewed as a problem of propagation along the optical 
axis, with periodicity corrections given by our newly obtained kernel (0). Section 5 
gives a more detailed explanation of the relation between discrete space and periodic 
structures, casting the propagator in terms of continuous variables flTOj) . We end with a 
brief summary. 
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2. Discrete propagator in one dimension 

Let us start with the free Schrodinger problem in discrete space and continuous time. 
We proceed by discretizing the derivative operator in the continuous wave equation. 
The corresponding wave function can be written as 4> n {T~) with n an integer denoting a 
point in discrete space and r a real parameter denoting time. Introducing the constants 
A (with units of energy x length 2 ) and a (lattice spacing), the Schrodinger dynamical 
problem is described by the equation 

- ^ [<Wr) - 20 n (r) + <j> n ^(r)} = irS^- (1) 
or, more concisely 

-\[^n + l(t) + ^ n - 1 (t)]=i^-, (2) 

which is obtained under the appropriate rescaling of time t = and the 

transformation ip n (t) = e lt (j) n (T). The eigenvalue problem corresponding to (T5]) can 
be solved analitycally by means of Bloch waves [8]: 

" o K+l + ^n-ll = E ^n, (3) 



^ = -Le* fcn , E k = -cosk, ke(0,2n) (4) 

This simple set of solutions gives rise to a spectral decomposition of the evolution 
operator 

u n ,ut) = r ' dwiwid**-*** (5) 

Jo 

and the solution of any initial-data problem satisfying fl2]) can be written as 



M*) = E u n>m (t)ip m (o). (6) 

meZ 

We are therefore in the position to write the propagator for positive times as K(n, m; t) = 
0(t)U n ,m(i), where it is left to compute the integral in (jSJ). We proceed by using the 
Jacobi- Anger expansion [9] of e~ lEkt leading to Bessel functions of the first kind, namely 

K(n,m;t) = —6(t) dke ik{n - m) e icoskt 



2tt w Jo 



—9(t) r dke ik{n - m) V i l e m Ji(t) 
2tt Jo lez 



6(t)t n - m J n . m (t), (7) 
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where we have used the Kronecker delta from the integration of the exponentials and 
the property J_ n = (— l) n J n in the last step. Noteworthy is the existence of a phase, 
which shall prove to be of utmost importance in the analysis of Feynman paths in 
further sections. The meaning of such a phase can be readily given by expressing ([7]) in 
polar form: K(n,m;t) = e i?rL//2 J„_ m (t), where L = n — m is a discrete length between 
endpoints. 

By one stroke, we can give the propagator for d- dimensional lattices of the cubic 
type: Let n, m be vectors denoting points in the lattice. The total propagator contains 
one- dimensional kernels appearing in the form of factors, leading to the overall result 

AT(n,m;t)=0(t) J[ /" "'•./„. ,„,(!). (8) 
j=i,...,d 

It is important to note that the propagators of more complicated lattices with nearest- 
neighbour couplings cannot be cast as the simple product above, a problem that shall 
be studied in a follow-up paper \±U\ for d = 2. In the following we verify the properties 
of K(n, m; t) by means of a few theorems on Bessel functions [9]. 

2.1. Properties of the new propagator 

The fundamental properties of continuous propagators can be verified in the discrete 
case as well. It is instructive to see that some of the well-known identities of Bessel 
functions have a meaning for our kernel. Recent studies on Bessel functions such as 
summation formulae pj] and operational methods [12] can be further utilized in the 
context of propagation. In this section we restrict ourselves to the most important 
properties. We finish by showing that the continuous limit of J7]) leads to the possibility 
of computing corrections due to a minimal spacing a. 

• Connection with the identity: It can be shown directly by evaluating the Bessel 
function using the ascending series. One has K(n, m; 0) = i n_m J n _ m (0) = 

•n—mx x 

b u n—m,0 u n,m- 

• Composition law: This shall be useful in a formulation of the propagator in terms 
of paths. The Neumann theorem for the addition of Bessel functions (see p. 358 
[9]) ensures this property. Let t x > 0,t 2 > 0. We obtain 

Y, K(n, I; h)K(l, m; t 2 ) = %" £ J^fa) Ji- m {t 2 ) 

= % n - m J n . m {t x + t 2 ) = K{n, m; h + t 2 ) (9) 

• The propagator as Green's function of the discrete Schrodinger equation: For t > 0, 
K(n, m; t) satisfies the homogeneous equation (j2J) by construction. To see this in 
terms of J n (z), we find that the recurrence relations for the derivatives of Bessel 
functions U' n [z) = J n _i(z) — J n+ x(z) (see formula (2), p. 17 in [9]) yield the result 
directly. When n = m and t is allowed to vanish, we recognize that 
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I dK(n n- 1) 

- [K{n + 1, n; t) + K(n - 1, n; t)} - i v ' ' ' = -iS(t)K(n, n; t) = -iS{t), (10) 

implying thus 

\ [K{n + 1, m; t) +K(n-l,m; t)} - = -iS(t)S n>m . (11) 

• Continuous limit: Here we restore our units by introducing again the quantity a. 
When a — > 0, we keep a(n — m) = constant = x — x' as our continuous variables, 
with the obvious implication that n — m must be large. In this limit, t ~ r/a 2 
becomes large as well, but as a consequence (n — m)/t = 0(a) — >■ 0. Therefore we 
invoke the asympotitc form of Bessel functions when both the argument and the 
index are large quantities (see C. VIII in [9J). This is the so-called approximation by 
tangents p3], where the argument of J is replaced by z/sec (3 such that f3 — > n/2. 
The corrections are under control by recognizing that the approximation can be 
expressed as the first term of a series through Meissel's second formula (p. 227 in 
0). We have 



(12) 

V27rz/(1 — sec z p) 1 / 4 
and using sec = z with z> 1, the Taylor expansions in 1/z lead to 



Jn—m {t) 



-iye lt exp 



2t 



V2ttU 



I ha 2 
2nirA 



exp 



,h(x — x') 2 

27A 



(13) 



Finally, we identify A = h 2 /m, where m is the mass of a quantum particle. The 
gauge factor e lt is due to our simplification in the discretization of the second 
derivative (the step leading to (J2J) from ([T])), whereas the factor of a is the 
infinitesimal length to be absorbed as a differential when we consider the summation 
of amplitudes, i.e. J2m a ^ / dx' . The factor of i m ~ n is cancelled by the prefactor 
in the last line of (GO). These steps lead to 



rw x r / m (.mix — x') 2 \ JN 

K{n, m; t) -* [ae>>] x ^— exp (>^-± J (14) 

which corresponds to the free propagator in continuous variables. It is amusing to 
see that the next order correction of Meissel's expansion leads to the first corrections 
of the propagator due to discreteness in space (or minimal length). This shall be 
investigated in the context of diffraction. 
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2.2. Connection with Feynman's sums over paths 

Now we investigate a formulation of discrete propagators leading to a summation over 
paths. The points to be touched upon are the existence of a phase associated to an 
action and a time- dependent weight function as in the usual Feynman path integral 

It is often pointed out that the original idea of composing propagators for small 
time intervals, comes from the fact that each propagator is a phase factor containing 
the classical action [HJ [151 EH]- In discrete problems, however, we do not have such a 
classical action and Feynman's recipe cannot be applied here straightforwardly. We 
recognize that the formalism of path integrals was designed to deal with problems 
in continuous space [17]. Nevertheless we may employ the composition property of 
discrete propagators ([9]) in order to investigate the process of multiple compositions 
along an arbitrary number of time slices. Denoting such a number as iV + 1 and defining 
vn+i = n,u = m,tN + i = t,t = 0, we can express our propagator as the sum of 
products 



N 

K(n, m;t) = E II K ( u j+u v v ~ t j) 

vi,...,u N j=0 
N 

= E IK" '-A',-: I'M,-- I,) (15) 
u\,...,u N j=0 

where the times are ordered as tj + \ > tj for each slice j. As a side remark, this 
expression is yet another representation of a Bessel function which, to the author's 
knowledge, has not been investigated in the mathematical literature. We proceed to do 
so, by taking equal time intervals tj+i — tj = t/(N + 1). We should note also that the 
kernel (J7J is symmetric K(n, m; t) = K(m, n; t) as a consequence of J- n (z) = (—) n J n (z) 
and of the inversion of the complex prefactor. Then, it is legitimate to replace 
i n ~ m J n _ m (t) = i\ n - m \ J| n _ m |(t) for each factor in the r.h.s. of (fT5l) . Furthermore, we 
express the complex prefactors i\ n ~ m \ in polar form to obtain 



K(n,m;t)= £ J[ e^~^ 2 J ]Vi+1 . Vj] {t/(N + 1)) 

vi,...,u N j=0 

(7T N \ N 

«•« e - v s \ n j\u J+l -u A it/{N + 1)) 
3=0 J j=0 

, N 

= Yl exp nV,(*/(^ + l)) (16) 

where in the last line we have identified SV+i,o = Y.f=o Wj+i ~ u j\ as the total length of a 
path joining the events (0,m) and (t,n). Similarly for Sj + ij. Therefore, the coherence 
of the sum is completely governed by the phase Sjv+i,o> which we are tempted to call 
action : The minimization of this functional gives a reduced class of paths containing, 
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in particular, the discrete straight line. One might call such a trajectory a classical 
path, although our dynamical problem has been defined exclusively for discrete quantum 
systems. The sum is dominated by the paths satisfying S m i n = n — m. We should note 
that the terms in the sum are also affected by a weight W = fL Js +1 - {t/N). Despite 
of being real, such a weight dictates the contribution of each term and is, at first sight, 
a path-dependent function. This is evident from the fact that each Bessel function 
contains explicitly the partial length Sj+ij in its index. In order to understand better 
this sum, we study the weights in the limit iV — > oo. Taking the first term in the 
ascending series of the Bessel function (formula (8) p. 40 in [9]), and denoting each 
path Vi,...,v n = {z/}, we can make explicit the separation between path-independent 
contributions w(t) and path-dependent contributions F[{z/}]: 

(. \ Siv+i,o N i 
WTT)) ^j^-W^Ml- (it, 

The final sum reads 



K(n,m;t)= ^ w(t)F[{u}] exp ( i-S N+1>0 ) 

Paths V 1 ' 

(18) 

and with this expression we find that Feynman's prescription applies to discrete problems 
with the appropriate amendments, i. e. a factor F correcting the weights. The partial 
path lengths are the key, as we have shown. The presence of F can also be analyzed 
on a conceptual basis: The paths contributing to our propagator are not in the same 
class as those appearing in continuous-variable path integrals. The latter turn out to be 
continuous curves, the majority of them being nowhere differentiable. In contrast, the 
paths that we are studying are discontinuous by assumption, jumping in Z (depicted 
in figured]). The paths that might approach dangerously to a non-measurable function 
in space-time are killed (or supressed) by the weight F. The paths that are favoured 
by F are those whose lengths do not approach to infinity and contain a finite number 
of discontinuities. And this is possible only if the continuous pieces of such curves 
are constant functions of t. This is a clear generalization of the integrals appearing in 
diffusive processes involving the Wiener measure [16], as well as a deviation from the 
traditional path integral formulation. It is, however, consistent with our mathematical 
procedures. 

To complete the circle, we finish this section by showing that taking (ITS]) as a 
starting point, we arrive again at a single Bessel function. First, we perform a change of 
summation index: Let {S} = Sjv+i.jV) S^o be the set of partial lengths characterizing 
a path {z/} of a fixed length S. We can separate the sum into 



oo 

E= E c(s)xj: 

{u} S=S min {S} 



(19) 
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Position 




Figure 1. Left panel: Three discontinuous paths of equal length joining (ti,Xi) — 
(tf,Xf). The blue path contains only one change in direction. Right panel: Comparison 
of a typical continuous path in Feynman's integrals (Red curve) and a discontinuous 
trajectory in our path formulation (Blue curve) 



where C(S) is a coefficient counting the number of paths of equal length. Obtaining 
C(S) is a combinatorial task and the procedure can be found in Appendix A. We have 



C{S) 



and noting that S„ 



S\ 



^ S+n—m ^ | ^ S+m—n ^ | 

= n — m, the sum (fTBl becomes 



(20) 



K(n,m;t)= ^ 



S\ 



S=n—m 



■ S+n-m \]( S+m-n \] eX P \^ 2 



7T 



s E 



{S} 



S at 



n 



(21) 



Moreover, we observe that the factor F is related to a multinomial coefficient, leading 
to 



E 



t 



s\ 



A' 



v §2(iv+i; 



{ 5 } :lcn gth =5V2(iV+l)/ njU(W 

Finally, the propagator can be written as an ascending series in t 



2, 



(22) 



K(n, m; t) 



, E ( 

S=n—m \ 

00 

jn—m ' 



cxp 



Of 5 ) 



s=0 

^ Jn—m\t) 



I)' 

S— ii-iii 1 ! , >'--/»-- n \\ ^ 2 / 

^ \ '2.1— n- m 

s + n — m)\s\ \2, 



t > 0. 



(23) 



We have shown that discrete quantum mechanics, in its simplest form (the free problem), 
allows a description in terms of paths and weights. The resulting paths are discontinuous 
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at a finite number of points in time and generalize the paths of diffusion processes. 
Moreover, the prescription we have obtained can be used again to recover a closed 
formula for the propagator — A book-keeping combinatorial procedure. 

3. Applications and new effects 

An important question about discrete and continuous quantum systems is related to the 
effects emerging in the presence of a minimal length. The obvious example in which we 
can distinguish clear differences comes in the form of diffraction, where it is known [7lfT8] 
that in the case of continuous variables, an infinite number of wavelengths is necessary 
in order to describe the propagation of edges (for example, diffraction in time). A 
direct consquence of space discreteness is the irrelevance of the Gibbs phenomenon in 
the Fourier decomposition of square packets. However we do not contempt ourselves 
with this simple conclusion: We extend further our study by describing the propagation 
of initial conditions such as the discrete version of the Moshinsky shutter [7] and the 
point-like initial condition. 

3.1. The expansion of a one- cell packet 

We start by analyzing the propagation of a packet of minimal width: one cell or pixel. 
In the continuous regime, such an initial condition cannot be conceived, but it is well- 
known that a packet with a square shape [19j EH] gives rise to non-trivial oscillations 
before the natural regime of expansion occurs. Let ip(0) be given by 



Contrary to the case of a point-like stimulus in continuous space, our evolved solution 
yields a wavefunction with non-trivial density in the position variable n: 



We can see the behaviour of the probability density in the plane (n, t) in figure [2j As 
the most prominent feature of this evolution, we find a regime of expansion which extends 
for all times. This is in sharp contrast with the continuous-variable problem, where 
infinitesimal oscillations typically appear around discontinuities in the initial conditions 
[T9~| [T8] . The edges are known to produce patterns approximately describable by the 
Cornu spiral, but the geometrical explanation provided by such a curve is not applicable 
to our discrete example. In fact, the regime of oscillations around an edge (the analog 
of Fresnel diffraction) is absent: We only find an expanding pattern with oscillations 
inside two wavefronts (bright diagonal fringes). The absence of the Fresnel regime is 
not the only notable effect. The natural frequency cut-off provided by the discreteness 



-0n (0) = 5 nfi . 



(24) 



\Mt)\ 2 = J 2 n(t). 



(25) 
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Figure 2. Left panel: Probability density in the plane n (abscissa) and t (ordinate) 
of a point-like source. We can see an expansion of the density at a constant velocity 
(set as unity). The expansion is accompanied by oscillations between the fronts n±t. 
These features cannot be found in the propagation of a point-like source in continuous 
variables. Right panel: Increased resolution, emulating a continuous variable n in the 
abscissa. 

of the problem (E max = A/a 2 ) suggests a maximal velocity of propagation, which can 
be taken as A /ha or unity in our conventions. The bright fringes of unit slope give rise 
to a conical structure in space and time, despite of the non-relativistic nature of the 
problem. The bright fringes indicate an abrupt change in the density, but not an actual 
discontinuity of the density function in the t variable — They appear due to a valley 
typically found in Bessel functions of large index. This peculiar slope or velocity shall 
appear again in the problem of the Moshinsky shutter. 

We note, in passing, that the aformentioned wavefronts can be related to caustics 
in the Bessel integral. These caustics can be easily determined by finding the stationary 
point with the phase function k(n — m) + tcosk. An appropriate description of the 
function around the caustics comes from the expansion of the Bessel function in terms 
of Airy functions [20] , covering the regions inside and outside of the conical structure 
(the Meissel expansions mentioned before cannot cover both regions simultaneously). 

3.2. The Moshinsky shutter: Diffraction in discrete space and continuous time 

Now we consider, for negative times, a free wave on the left semi-axis moving to the right. 
A blocking (absorptive) screen at the position n = prevents the wave to propagate to 
the positive axis. In order to mimic Moshinsky's problem, such a wave must be taken 
as a Bloch wave, which solves the problem without boundaries. At t = the shutter is 
removed and the evolution studied: 
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ikn 



^ n (0) = 0{-n)e 



m=—oo 
oo 

i ^ ] Jn+m{t)& 
m=0 

oo 



5(0) = 1 

„(t)e* fcm 

—i(k—n/2)m 



(26) 



(27) 



— E Jn+m{t)^ m 
m=0 

Although the series above has the form of a generating function with complex parameter 
£, there is no closed form reported in the literature for this expansion (the summation 
index does not take negative values). However, we can extract some properties of this 
solution by exploiting some identities. For example, we can find the value of if) at 
position n when the wave function at the shutter position is known, i.e. 



ikn 



n-1 



-i(k— n/2)m 



m=0 



(28) 



where the limit of the sum becomes finite. The function ipo(t) has notable properties, 
for instance 



Mt)\k= g -n/2 + Mt)\k=./2- q = e ltsinq + J (t). (29) 
A quick proof of this formula can be obtained by direct substitution of the second line 
of ( 1271) on the l.h.s. of (|29l) . Then by means of the Jacobi- Anger expansion we have 

Rising = + £oo =q jjfi ( e mg + (_)n e -tn«) and the r h g of Qfify f H ows . 

A particular case of fl29l) (q = vr/2) gives the possibility of extracting ip for k = 
in closed form: 



V>n(*)| 



fc=0 



~e u + J (t) 



n-1 



(30) 



m=0 



Another limit of interest is the large time behaviour of if). This can be estimated in terms 
of the Meissel asymptotic form of Bessel functions as long as (n, t) is not close to the 
caustic. We write J n+m (t) w -y/^exp (P n+m + iQ n+m ) where Q and P are well-known 
functions (formulae (1) and (2) in p. 228 of [9J). Upon substitution in ( 1271) 



/ 2 00 

if) n (t) « A/ — E exp (P n+m + iQn+m) ■ 

V m =0 

leading to the typical decay of the density as the inverse of t 



(31) 



TXt 



exp (P n+m + iQ 



n+m J 



(32) 



m=0 



For a description of the probability density without approximations, see figure [3J 



Exact propagators on the lattice with applications to diffractive effects 



12 




Figure 3. Probability density for the problem of diffraction in discrete space 
(abscissa) and continuous time (ordinate). In the first panel we show the position 
of the shutter denoted by S, and a pair of lines indicating the propagation of fronts 
at approximately constant velocity; they propagate both forward and backwards. The 
value of the Bloch quasi-momentum k is varied for each panel and it is indicated on 
the top of each plot. The six shown patterns repeat in reversed order as k increases 
beyond 7r/2. These effects can only happen in discrete space: An ever increasing k 
does not lead to faster propagation. Notably, the structure outside the cone t±2n/e 
remains intact, as opposed to the richier structures found in the continuous variable 
solutions (Moshinsky functions). For this computation we have used the expansions in 
(|27| up to 30 terms. 



3.2.1. Maximal velocity in a non-relativistic framework This interesting property is a 
reminiscent of the behaviour seen in the point-like case. It can be derived fairly easily 
by considering short times and arbitrarily large index, such that each Bessel function in 
the expansion ( 1271) can be replaced by the first term of the ascending series: 



m=0 



t 



n+m 



(n + m)\ 



+ ■ ■ 



-i(k— n/2)m 



(33) 



Keeping lowest order in powers of t, the probability density can be approximated by 
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^•"(SfG)"- (34) 

we can use the Stirling approximation to obtain constant probability density curves (say 
\%l)\ 2 = a) in the form 

/2a 1 /(2n)\ 

t = ±i — - — In. (35) 

The cones described by these curves (straight lines for a = 1) are shown in the first 
panel of figure [3j The presence of the cones is independent of the quasi-momentum 
k. This property is found exclusively in discrete problems, as non-relativistic quantum 
mechanical systems in continuous variables do not suffer the restrictions of a maximal 
velocity of propagation. 

3.2.2. Circular property of the solution The time-diffractive solution (|28|) can be 
related to circular waves. This property can be recognized by inspecting the expansion 
( 1271) . When we replace the time and quasi-momentum variables by radius and angle in 
the form 1 1- > r, k i— > (p, we immediately see that 



^(r) = e m ^$ m+n (r )V ) (36) 

m=0 

where $j(r, ip) is a circular wave of angular momentum / G Z and unit energy E = 1. 
This means that the function ty n (r,ip) = e~ inip if} n (r) satisfies the Helmholtz equation of 
unit wave vector in polar coordinates x = t cos k,y — t sin k, namely 



V| )B * n (r,^)+* n (r,^) = 0, \/n (37) 

and it corresponds to a free wave in two-dimensional space with a non-vanishing 
boundary condition at the origin: 



,1 for I < 

¥i(°.¥0 = < n +u . (38) 
otherwise 

i.e. a left-handed solution, since positive angular momenta do not contribute in (138]) . A 
depiction of the diffractive process in terms of circular wavefronts is shown in figure HI 
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Figure 4. Diffraction in discrete space as a circuiar wave. Each panei represents 
the probability at a fixed site, going from dark areas (low probability) to bright areas 
(high probability). The radius stands for t = \J x 2 + y 2 , while the angle denotes the 
quasi- momentum k = arctan(j//x). On top of each plot we indicate the value of the 
position n. A circular pattern around the origin emerges, showing that a signal arrives 
at the specified point after a fixed time, irrespectively of the initil momentum of the 
wave. One of such circles is depicted in the central panel of the lower row. 



3.2.3. Recursion in the position variable and algebraic properties Since our time- 
diffractive solution is related to circular waves, we expect a number of well-known 
properties to be applicable in our study, in particular when we try to evaluate the 
wave function at a certain position by knowing its value at a distant point. We can 
extract a recursion relation in n by noting that a circular wave of arbitrary angular 
momentum can be obtained by "cranking up" the solutions of lower angular momenta. 
In the language of Bessel functions, this corresponds to a ladder operator which is linear 
in the radial derivative. For the full circular wave, this process can be generated by a 
chiral operator. We have 



-cL ± d r 
r 



$„(r,^) (39) 



leading to the validity of ( )39|) for any superposition of $'s, in particular 
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$f n+1 (r,(p) = - e 



up 



(40) 



Its application to \l/o an arbitrary number of times gives 



*o(r,y»). 



(41) 



Finally, using the definition of \& in terms of ip, we obtain a new representation of the 
wavefunction at distant points in terms of the solution at the shutter: 



il> n [t)\k 



Akn 



-d k + d t 



Mt)\i 



(42) 



With this, we establish that "hopping" over discrete sites provides means to evaluate 
diffractive solutions. This process is related to the algebraic properties of the problem, 
obtained from the factorization of V^ fc . 

3. 3. An edge producing a Bessel profile 

Many other examples of evolution can be studied in exact form by using several 
results on the closed summation of Bessel functions. In connection with diffraction 
problems by edges, we envisage a wide class of results by means of degenerate addition 
theorems where the summation index starts at 0. Lommel expansions constitute another 
possibility. Take, for instance, an initial condition in the form of a ramp : 



^n(O) = 9{n) 



111 



(43) 



This poissonian profile represents the wave function of a particle that is trapped near 
the shutter, with a fastly decreasing tail. The function vanishes abruptly for negative 
positions and it can be regarded as a member of the family of edges. The real parameter 
b determines the position of the maximum, e.g. for 6=1 the maximum is exactly at 
the position of the edge. The evolution is obtained by computing the sum 



1 Jn—m\t) 



T — 

[i n {t - 2ib) n/2 



t- n/2 J n (Vt 2 - 2ibt) 



(44) 



where the second factor is manifestly regular at t = 0. In this derivation we have used 
the Lommel expansion (formula (5) p. 141 



j l/ (zVT+^) = (i + wy/ 2 J2' 

m=0 



-J, 



ml 



v+rn \ 



(45) 
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From this result we can infer many properties that coincide with our previously obtained 
effects. For instance, using the first term in the ascending series of J one can find 
constant probability curves involving n and t, finding a short time behaviour of the 
caustics and a maximal velocity effect. A decreasing density in the form 1/t can be 
found from the asymptotic form of J 2 . When 6 = 0, we recover the point-like source. 



4. Scattering by an edge in discrete space as a problem of propagation 

It is well-known that the shutter problem in continuous variables represents a diffraction 
problem in two spatial dimensions — Fresnel diffraction. Is it possible to show a similar 
analogy between our discrete problem and a scattering system in two dimensions? In 
the rest of this paper we clarify the connection of our evolving edges with a genuine two- 
dimensional problem described by the Helmholtz equation. To this end, we demonstrate 
how a problem of scattering with a periodic structure as a background (one discrete 
variable and one continuous variable) can be written in terms of a propagator. The 
corresponding kernels can be given in explicit form. 



4-1. Propagator for edge diffraction in two continuous variables 

For the sake of completeness, let us deal first with the two-dimensional continuous 
problem of a scalar plane wave scattered by an edge. The problem can be written in 
terms of forward and backward propagation (transmission and reflection) in the following 
manner: Suppose that a plane wave arrives at a screen positioned at z = and parallel 
to the x axis (at this point, we make no further assumptions on the nature of such a 
physical obstacle). We have ip(z < 0, x) — exp(ik z z+ik x x) with k z , k x the corresponding 
components of the momentum. At z = we can impose an arbitrary boundary condition 
determined by the form of the screen; for definiteness we may put ip(0, x) = 6(— x)e x . 
Set the energy E = k 2 + k 2 such that 



V 2 + E 



i>(z,x) 



0. 



(46) 



Since the solution of the problem satisfies the Helmholtz equation, it can be written 
formally as a continuous superposition of the form 



ip(z, x) = A R / dk(p^ exp (ikx + izV E — k 2 ) + A L / dk(j)\ exp (ikx — izV E — k 2 J 

J — CO ^ ' J — OO ^ ' 

(47) 

where L and R indicate solutions propagating to the left and to the right, respectively. 
Now, let us assume no reflection (A L = 0, A R = 1) 



dk<j)k exp (ikx + izV E — k 2 ^j 



(48) 
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and let us evaluate ip at z = in order to determine the distribution in k: 



V(0,ar) 



dk(p k e 



ikx 



dkip(0, x)e 



-ikx 



1 '' x 

Z7T J-oo 

Substituting the inverse Fourier transform of the initial condition we have 



(49) 



tp(z,x) = — / dx'i/j(0,x') / dkexp Uk(x — x ) + izy E — k 2 

ZlX J—oo \_J — oo ^ 



(50) 



Finally, the propagator in fictitious time z and parameter i? is given by the factor in 
brackets 



\ roo . . 

K E (x, x'; z) = — / dk exp (ik(x — x') + izV E — k 2 ) . 



(51) 



This integral can be evaluated upon a change of variables [21 J and the appropriate choice 
of an integration contour. Let rj = E~ l / 2 {k+i\/ E — k 2 ) and r? -1 = E~ 1 / 2 (k—i\ / E — k 2 ), 
with an integration line given by the segments {Re(^) < —1} U {77 = e tS \0 < 5 < tt\ U 
{1 > Re(?7) > 0}. As there are no poles in the way, the upper semi-circular curve and 
the segment (0, 1) can be continuously deformed such that r] G (— 00, 0). The integral 
in floTT) is transformed to 



K F (x,x': z) 



VE r° 



4tt 



drj 1 — 7] 



-2 



exp 



i](x — x + %z) + 7] 1 (x — x — iz) 



drj <r] 



47T JO 

The closed result is 



-2 ^^■^r)(x~x'+iz)+ri- 1 (x-x'-iz)^ _ ^(x-x 1 +iz)+r]~ 1 (x-x'-iz)} 



K E (x,x'] z) 



-izVE 



2^J(x - x') 2 + z 2 



(i) 



EJ(x - x') 2 + z< 



(52) 



(53) 



where we have used the integral representation of the Hankel functions H±l [13] and 
trivial algebraic manipulations. The result found here should not be confused with the 
usual Green's function of the 2D Helmholtz equation subject to Sommerfeld's outgoing 
radiation condition (given by the Hankel function Hq 1 ^ alone). Our boundary condition 

has been defined along the x axis at z = 0. 

It is important to note that in fl53|) . the two limits J{x — x') 2 + z 2 ^> l/yl (short 



wavelengths) and (x—x') 2 /z 2 1 (paraxiality) must betaken simultaneously in order to 
recover the free Schrodinger propagator by means of the asymptotic form of the Hankel 
function: 
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The prefactor in brackets is related to the planewave propagating originally to the right 
on the left semi-plane z < 0, with a wavelength A = 2n/\/E. 

Although the problem of diffraction by an edge has been discussed exhaustively in 
the literature, no description of the propagator (|53|) can be found in previous works. We 
should stress that this is the Euclidean version of a relativistic problem described by the 
Klein-Gordon equation [2TJ. Since the solutions to relativistic diffraction in time were 
identified as Lommel functions [7], we should expect similar results in our Euclidean 
case. 

In connection with the validity of the approximations leading to the Schrodinger 
kernel, it is worth to mention that the problem of diffraction in propagation axis z 
is characterized by infinitesimal oscillations near z = due to a dependence of the 
form (x — x') 2 /z in ( |54l) . On the other hand, the kernel (l53l) exhibits a dependence on 
\J(x — x') + z 2 which is regular at z = and does not produce infinitesimal oscillations. 
Now that we have set the continuous example, we move to the discrete problem. 

4-2. Continuous propagation axis and discrete transverse dimension 

We turn our attention to the Helmholtz equation of one continuous variable z and 
one discrete variable n. This problem also corresponds to a wave propagating in a 2D 
medium with one periodic direction (see section |5] and figure E]). The discrete derivative 
operator along the variable n contains twice the identity operator, therefore it can be 
simplified by a redefinition of the energy in the full Helmholtz equation, leading to the 
analog of the l.h.s. of (J2J). We have 

+ Vwi(*) + ^n-i(z) + Eij n (z) = 0. (55) 

In complete analogy with the continuous 2D example, we state our boundary 
problem as follows: A Bloch wave propagates in the left semi-plane towards an 
absorptive edge, i.e.ip n (z < 0) = e lk ^ n + tk ^ z ^ where k x is Bloch's quasi-momentum and 
k z is the momentum corresponding to coordinate z. The energy E = k 2 z — 2 cos k x is 
fixed. The most general superposition which solves (1331) when the wall is absorptive can 
be written as 



ip n (z) = J dkfaexp (ikn + izy/ E + 2 cos k^j , z > (56) 

where the integration over k extends to the first Brillouin zone (0,2tt). The boundary 
condition at the position of the screen implies 



r 2n 1 

^ n (0) = / dkfae ikn ^<p k = —Y, e- ifc "> m (0). (57) 

JO 2,71 rj 
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n=3 




n=2 




n=l 




n— u 




z 


Medium 1 




Medium 2 


^^^Z Screen 











Medium 1 

/ 




n=0 



n=l 



n=2 



Figure 5. Left panel: A periodic background realized through the alternation of 
two materials (solid state) or two potentials (quantum case). A screen blocks a wave 
propagating along z. Right panel: The two media represented as potential barriers 
along the x variable. A sketch of the profile for localized states is shown, with coloured 
areas indicating the overlaps and nearest-neighbour interactions. 



For definiteness, we may choose 



^n(O) 



ik x n 



9(-n), 



-I CO 

__ p i(k-k x )m 
27T Z ^ 



J(k x -k)/2 



(5f 



m =o 2nism[(k x - k)/2] 

although the propagator we are seeking for does not require such an explicit function. 
Using the Fourier series in (|57|) for we have 



— V Vm(0) [ 1^ dke ik{n - m)+izVE+2cosk 



(59) 



and, as before, the propagator can be read-off from the expression between brackets: 



K E (n,m;z) = — 

Z7T 



f 2n ' 

/ dke^ m)+iz\/ E+2 cos k 

Jo 



(60) 



The discrete kernel ( |60l) can be expressed in compact form by employing a trick used by 
Dattoli [T2] (with reminiscence from a formula by Filon, see formula (5) p. 51 [9J). A 
translation operator of variable E and parameter 2 cos A;, acting on the function e tz ^ 
can be expanded in terms of Bessel functions as follows: 



e izV E+2 cosk = exp 



zcos*(-2zA 



(61) 



where each Bessel operator is understood in terms of its ascending series. Inserting back 
in ( 16 Op and performing the integration over k one finally gets the simple expression 
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K E (n,m; z) 







i Jn—m ( 2? „ „ I 6 



izVE 



dE 



(62) 



which ressembles our discrete Schrodinger propagator ([71) in the form 



K E (n,m; z) = K Schl - 6dingcr In, m; ~ 2i -^j e ' 



zVe 



(63) 



This operational form has several uses, according to our needs [22] • For example, 
consider the following identity for some analytic / 



z^E _ e izVE 



[1] 



(64) 



where 1 is the constant function equal to unity In Appendix B we show that this 
identity can be approximated (in a short wavelength regime) by the expression 



/ 



/ 



_|_ %Z j" 



(65) 



This allows to compute the first corrections to our propagator Kg by replacing the 
discrete Schrodinger kernel K in / and /": 



Ke(ti, m\ z) w e ! 



K I n,m: 



iz 



16E 3 / 2 



K [ n + 2,m: 



E. 



+ K [ n — 2,m: 



+ 2fT n, m; 



(66) 



'£/ V VEJ \ VEJ 

The first term in (166"1) establishes the correct limit for short wavelengths. It gives the 
discrete Schrodinger theory with fictitious time t = zj y/~E and the additional plane wave 
factor e lzyJ ^ '. For diffraction problems, we can establish a similar approximation for the 
wave function in terms of the Schrodinger solution ip in (1281) 



z\fE 



IZ 



E, 



x W 



n+2 



16E 3 / 2 
+ 1pn-2 



X 



2^n 



(67) 



'Ej WE J WE J 

This is our final result: The wave function of a diffractive experiment in 2D ressembles, 
in a suitable approximation, the problem of diffraction in discrete space and continuous 
time. The parallelism Shutter ■H-Fresnel between the continuous problems is also valid 
for our discrete settings. Moreover, the effect is accompanied by the short wavelength 
corrections analogous to those separating the Hankel function (1531) from the gaussian 
kernel (1341) . This has consequences in the propagation of waves on a periodic background 
(for example, a photonic structure). 
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5. Connection with tight-binding arrays and photonic structures 

The connection between the discrete version of the Schrddinger equation and tigh- 
binding arrays is rather direct. Let us focus in the lowest energy band of a homogeneous 
polymer chain with nearest neighbour interactions. In the appropriate units, the 
hamiltonian for this band is 

H = ~\ E [\n)(n + l\ + \n + l)(n\] (68) 
neZ 

where any localized wavefunction £ at site n can be obtained via the kets above: 
£ n (x) = (x\n). Our present work on the propagator leads to the solution of the time 
dependent-problem described by the wave function 

ij>(x,t) = £ (69) 

neZ 

with the continuous-variable propagator restricted to one band given by 

K poXymcr (x,x'\t) = ^n(x)^(x')K(n,m;t). (70) 

n,msZ 

In a polymeric problem, we are usually given the localized functions £ n (a?) describing a 
quantum particle at isolated sites. When we deal with the propagation of some initially 
known packet ipo(x), we assume that the packet can be completely resolved by the basis 
of localized states ipo(x) = J2 n C n ^n(x) (Wannier). In other words, we rely completely 
on the information given by the overlap / dxifjo(x)^(x). This is the first step towards a 
time-dependent description of conductive (conjugate) polymers of a more complicated 
nature [2]. 

6. Summary 

In this paper we have calculated a propagator in discrete variables, with no precedent 
in the standard literature jTU [151 US]. In order to understand such a novel object, we 
have studied its properties and extended the Feynman path sums to discrete variables. 
We discussed some relevant examples, including the diffraction by edges and the effects 
emerging from a minimal spacing. We established the mathematical form of the solutions 
and gave a detailed comparison with a problem in two dimensional space in a periodic 
background. A possible realization has been proposed in tight-binding arrays. The wide 
interest in photonic structures suggests applications of our results in this well 
as solid state physics in time domain. 
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Appendix A. Paths of a fixed length 

Here we find the number of paths of a constant lenght S. We are looking for all 
the solutions of the equation Sjv+i,o — S — constant, subject to the constraints 
uq = m, vn+i — n - These conditions lead to configurations in which the segments of the 
paths can change directions, displacing from a point to its neighbours discontinuously. 
For each path we can define Si as the number of total steps going left and S r for steps 
taken to the right. We obviously have 



S=Sl + S r (A.l) 

Such steps can be taken in any order, but conspire to give the same total length 
S = Sjv+i,o and the same initial and final position. Counting them is equivalent to find 
the number of combinations of S r segments in S places (equivalently Si in S places). 
As the partitions matter but not the permutations among fixed segments, we have the 
binomial distribution 



Now, S r turns out to be a fixed number, independent of the number of time slices N. It 
corresponds to the maximum displacement to the right in a path with a single change 
in direction (see the blue path in the left panel of figured]): S r = h(S + n — m) which is 
always an integer, as S and n — m must have the same parity. From this condition we 
also obtain Si. Therefore we have 



C(S) = 




(A.3) 



With this calculation, we can find the Feynman kernel as indicated in section 2.2. 



Appendix B. Derivation of an operational formula 

Here we derive the approximation (|65|) . Consider the general identity 



for / and g analytic and [1] the constant function. This identity can be derived easily 
by noting that 
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^ = 9(x) ( j- + [1] (A. 5) 

dx w \dx g(x) J 1 1 K ' 

and similarly for its powers. For simplicity define now h(x) = g'(x)/g(x). Taylor- 
expanding / on the r.h.s. of ( ]A.4j) leads to 




(A.6) 



Now we expand the binomial (d/dx + h(x)) n and keep the only term without derivative 
operators (they act on the constant function [1], therefore they vanish). The resulting 
expression contains terms of the form h(x) a d^.h(x) c with a + b + c = n — 1. From these 
terms we only keep first derivatives, i.e. b = 1. This corresponds to small inverse powers 
of the energy E in ( )65l) whenever x = izyfE and g(x) = e x . This approximation leaves 
us with the terms 



"~ 2 n(n — 1) 

£ h(x) j d x h(x) n ~ j = — — -h(x) n - 2 ti(x) (A.7) 

3=0 2 

as the only contribution for each term in the Taylor expansion. Substituting this back 
in (1A.6|) yields 



E /n( .° )w(n 9 1] h(xT-*h'(x) = ^P-f'(h(x)). (A.8) 
Finally, we have 

With this expansion we can give certain limits for the operational form of our propagator 
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